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Flux change associated to metabolite production correlates with genetic rewiring 
in catabolic but not in biosynthetic metabolites 

We summed the metabolic fluxes contributing to the production of a particular metabolite (m) 
in WT conditions (£™ t ), and in each of the backgrounds (bgi) defined by deletions of single 
genes (S^). We divided both measures (each one normalized by their corresponding flux of 
the biomass reaction, that we denote H^j. and S^M to compute the fold-change of metabolite 
m production in all backgrounds (i.e., Hj^/E^). Fold-change in most metabolites significantly 
correlates with genetic rewiring (table SI, this is consistent with the association between genetic 
rewiring and metabolic readjustment; see also main text). 

Notably, correlations are stronger for redox equivalents and currency metabolites (nadh, 
nadph, h, phosphate, atp), the serine —> glycine —> formate pathway intermediates (e.g., thf, 
mlthf, for, phosphoserine, serine, glycine), oxidative phosphorylation metabolites (e.g., quinones), 
glycolytic intermediates (e.g., glucose, glucose-6-phosphate, f6p, g3p, pyruvate) and TCA cycle 
intermediates (alpha ketoglutarate, isocitrate, oxaloacetate) . Conversely, most intermediates 
from diverse biosynthetic pathways either correlate more weakly or do not correlate at all with 
genetic rewiring. 

Catabolism in the Sacharomyces cerevisiae iND750 model 

Catabolic action consists of the breakdown of nutrients into energy (usually in the form of 
ATP), reduced equivalents (e.g., NADH, NADPH), and smaller components to be involved in 
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biosynthetic reactions (precursors). All of this can be produced in the oxidative pathway af- 
ter initial glucose uptake, through glycolysis and citric acid cycle (or tricarboxilic acid cycle, 
TCA). In addition, the most efficient ATP production pathway in aerobic conditions is oxidative 
phosphorylation that is able to couple the oxidation of reduced equivalents (produced in gly- 
colysis and TCA cycle) to ATP synthesis. Thus, a considerable proportion of the flux through 
both glycolysis and TCA cycle is devoted to the supply of reduced equivalents to oxidative 
phosphorylation. 

Beyond these "canonical" catabolic pathways, the iND750 model also predicts artifactu- 
ally the production of an important fraction of ATP. Specifically, part of the glycolytic 3- 
phosphoglycerate is diverted towards serine and then glycine that later becomes fermented to 
C02 by the glycine cleavage system. This cycle also produces formate that is then oxidized to 
C02 while at the same time reducing quinone, thereby also feeding oxidative phosphorylation. 
The model predicts this pathway as one of the main catabolic ones. For instance, this pathway 
takes over ATP production (accounting for about 60% of the total; see table S2) when oxidative 
phosphorylation is compromised (e.g., in the ACOX1 background). Note as well that the path- 
way includes genes from several annotation groups that would not be expected to contribute 
to glucose catabolism, such as the glycine/serine (e.g., SER1, SER2), or the folate group (e.g., 
GCV3). Certainly, this form of "glycine fermentation" is artifactual in yeast, and we found only 
one reference to a similar metabolic process occuring in Clostridium purinolyticum [1]. This 
catabolic role determines the rewiring and stability properties of the associated genes, and in- 
deed they act in a similar way to the rest of "canonical" members of catabolic pathways (e.g., 
there is an abundance of genetic interactions between SER1 or SER2 with glycolytic, TCA cycle 
or oxidative phosphorylation genes; see Fig. 3, main text). 

Other genes that belong to biosynthetic annotation groups also displayed epistatic interac- 
tions characteristic of catabolism (i.e., weak interactions with other genes of the catabolic core), 
such as AAT1, ALD6, or ALT2 (all associated to the Tyr, Trp and Phe group). This is due to 
their role in inter-conversions important for catabolic balances, such as a-ketoglutarate f± ox- 
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aloacetate, acetate f± acetaldehide and a-ketoglutarate f± pyruvate, respectively. On the other 
hand, they also exhibit several synthetic lethal interactions that are linked to their biosynthetic 
roles (more details below). 

Genetic interaction stability as a function of sign, strength and functional 
association 

In Fig. SI, we show the difference in interaction stability for different types of interactions. 
Interactions between genes in the same metabolic annotation module are more stable than 
those between different annotation groups; strong interactions are more stable than weak, and 
negative interactions are more stable than positive. To eliminate possible confounding effects 
between these factors, we performed separated Wilcoxon tests controlling for each of these 
factors, generally confirming the independent influence of all three on stability (table S4). 

Case study: Metabolic reorganization and genetic rewiring in the AMIR1 
genetic background 

In order to better understand the metabolic underpinnings of genetic rewiring, we analyzed in 
detail the changes occuring in the MIR1 background, the mitochondrial phosphate carrier needed 
to provide phosphate for ATP synthesis through phosphorylative oxidation. This gene represents 
a typical example of genetic hub, with 31 epistatic interactions (24 of these being weak) that 
encompas 16 different metabolic modules (many of them are clearly biosynthetic, e.g., Arg/Pro 
or Thr/Lys). 

In Fig. S3, we schematically depict some of the relevant metabolic changes that occur after 
deletion of this gene. The impairment in phosphate import to mitochondria caused by the 
deletion is "solved" in first instance by the antiport of phosphate with malate by the gene DIC1. 
Malate is then transformed to oxaloacetate in order to be reimported to mitochondria. This 
process results in a net extrusion of redox equivalents to cytosol that works as a kind of "redox 
shuttle". Ultimately, many metabolic pathways readjust to buffer the generated imbalance 
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(including some biosynthetic ones that conceal redox-balancing possibilities). Specific examples 
include: 

• Proline starts being synthesized from orotate in the cytosol to be then imported to mito- 
chondria, where it becomes tranformed again into orotate through several steps. The entire 
process transports redox equivalents from cytosol (NADPH) to mitochondria (NADH) and 
converts an additional molecule of mitochondrial NADPH into NADH. 

• The glycine cleavage system (artifactually) stops synthesizing ATP and redox equivalents, 
avoiding further release of NADH and NADPH to the cytosol. As a consequence, respira- 
tory activity increases as compared to the WT. 

• First steps of lipid synthesis from DHAP switch to an alternative pathway using NADH 
instead of NADPH. 

• Further reorganizations (Fig. S3). 

How are all these changes reflected on the genetic landscape? Firstly, activated compensatory 
mechanisms (e.g., mitochondrial malate synthesis and its cytosolic conversion to oxalacetate) 
start interacting positively with those genes experiencing the buffering changes (i.e., genes in 
the oxidative phosphorylation apparatus) but also among them, as part of the new optimal 
solution (Fig. S4). Second, since different compensatory mechanisms often contribute to similar 
functions (e.g., importing redox power to mitochondria) they can in turn interact negatively 
interacting among them (e.g., PUT1 and ALD6). Finally, newly revealed negative interactions 
can reflect as well third-instance alternatives able to compensate for the loss of some genes in 
the new solution. These types of genetic rewiring are summarized in Fig. S5. 
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Table SI. Fold-change of metabolic production significantly correlates with rewiring, a signal 
that appears particularly strong in catabolic metabolites (}P < 10~ 20 , $P < 10~ 10 , *P < 0.01). 
Metabolites are colored according to broad functional context where they appear. 
Red-catabolism; I ilue- biosynthesis; violet— currency metabolites and related molecules (e.g., 
H + or phosphate); yellow-metabolites belonging to artifactually catabolic pathways (see 
below), v -other types fo metabolites (such as water, inorganic nutrients, etc.). 



ATP synthesis mechanisms 
ATP synthase FTHFL PGK + PK 

WT 100 (60%) 26 (16%) 37 (23%) 

Background 

ACOX1 0 54 ( 60%) 36 (40%) 

Table S2. ATP production (mmol grDW^ 1 h^ 1 ) by different mechanisms in WT and in 
ACOX1 background (percentage of the total ATP produced in the cell in each case in 
parenthesis). FTHFL: formate tetrahydrofolate ligase, PGK: phosphoglycerate kinase, PK: 
pyruvate kinase. See second section Supplement for details. 
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Table S3. Distribution of genetic interactions in catabolic modules. 
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ALT2 


MET22 


0.5 


AlaAsp 




Cys 


ALT2 


FUM1 


0.5 


AlaAsp 




OxPhos 


GLY1 


MIR1 


0.5 


ThrLys 




MitTrnsp 


LPD1 


PPK1 


0.5 


TCAcyclc, GlySer 


, Glycls 


Glycls 


LPD1 


RPE1 


0.5 


TCAcyclc, GlySer 


, Glycls 


PentPh 


LPD1 


FBA1 


0.5 


TCAcyclc, GlySer 


, Glycls 


Glycls 


SER2 


PFK1 


0.5 


GlySer 




Glycls 


SER2 


RPE1 


0.5 


GlySer 




PentPh 


SER2 


FBA1 


0.5 


GlySer 




Glycls 


SER2 


AAT1 


0.5 


GlySer 




AlaAsp, ArgPro,TyrTrpPhe 



Table S3, (continued) 
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PFK1 


SERl 


0.5 


Glycls 


pl..O TD ,J,. 

blyocr.rax 






0.5 


rcntrli 


MitTrnsp 




SER1 


0.5 


rcntrli 


pi„Q„ r TD ,J „ 

ialyocr,rax 


MIR1 


A ATI 


0.5 


MitTrnsp 


AlaAsp, ArgPro,TyrTr pP he 


FBAl 


SERl 


0.5 


Glycls 


lalybcr.rdx 


A ATI 


SERl 


0.5 


AlaAsp, ArgPro, TyrTrpPhc 


ni..C TD J,. 

Ulyocr,rax 






0.5 


AlaAsp, ArgPro, TyrTrpPhc 


TyrTrpPhc 


COX1 


ATP 8 


1 


Uxrnos 


Uxrnos 


COX1 


COB 


1 


OxPhos 


OxPhos 


COX1 


PGI1 


1 


OxPhos 


Glycls 


COX1 


IDP1 


1 


OxPhos 


'IV"' A ...... 1 . . 

1 UAcyclc 


COX1 






UXr IIOS 


MitTrnsp 


ATP 8 


COB 


1 


OxPhos 


OxPhos 




GCV3 




Uxrnos 


I'olatc.Ulyocr 


ATP 8 


PGI1 


1 


OxPhos 


Glycls 








Uxrrios 


1 U-Acyclc 








uxrnos 


GlySer 








Uxrnos 


MitTrnsp 








Uxrlios 


AlaAsp, ArgPro, TyrTrpPhc 








Uxr DOS 


vjlyocr , Pclx 








Uxrnos 


TyrTrpPhc 


COB 






Uxrnos 


Glycls 


COB 






Uxrlios 


TCAcyclc 


COB 


MIR1 


1 


OxPhos 


MitTrnsp 


GCV3 


LPD1 


1 


Folate,GlySer 


TCAcyclc, GlyScr, Glycls 


PGI1 


LPD1 


1 


Glycls 


TCAcyclc, GlyScr, Glycls 


IDP1 


MIR1 


1 


TCAcyclc 


MitTrnsp 


IDP1 


ALD6 


1 


TCAcyclc 


TyrTrpPhe 


SER2 


MIR1 


1 


GlyScr 


MitTrnsp 


SER2 


SERl 


1 


GlySer 


GlyScr, Pdx 


PFK1 


FBAl 


1 


Glycls 


Glycls 


MIR1 


SERl 


1 


MitTrnsp 


GlyScr, Pdx 


MIR1 


ALD6 


1 


MitTrnsp 


TyrTrpPhe 



Table S3, (continued) 



Intra-group 2.7 2.3 X 1CT 5 

Strong 

Intcr-group 5.1 



Weak 


Intra-group 
Inter-group 


8.3 
8.7 


0.97 


Strong 


Positive 
Negative 


7.1 
2.4 


2.1 x 10~ 9 


Weak 


Positive 
Negative 


10.9 
7.1 


1.8 x 10~ 7 


Intra-group 


Positive 
Negative 


7.1 

4.2 


0.013 


Intcr-group 


Positive 
Negative 


10.3 
6.0 


1.78 x 10- 1 " 


Intra-group 


Strong 
Weak 


2.7 
8.3 


9.7 x 10~ 8 


Inter-group 


Strong 
Weak 


5.1 
8.7 


1.9 x 10-7 


Positive 


Intra-group 
Inter-group 


7.1 
10.3 


0.014 


Negative 


Intra-group 
Inter-group 


4.2 
6.0 


0.002 


Positive 


Strong 
Weak 


7.1 
10.9 


1.8 x 10~ 4 


Negative 


Strong 
Weak 


2.4 
7.1 


5.2 x 10~ 18 



Table S4. Genetic interaction stability as a function of sign, strength and functional 
association. 
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Gene 


Pleiotropy 


Annotation 


TPIl 


41 


Glycls 


PFK1 


41 


Glycls 


FBA1 


41 


Glycls 


C0X1 


39 


OxPhos 


ATP8 


39 


OxPhos 


COB 


39 


OxPhos 


MIR1 


34 


Mit . Trnsp 


PGI1 


30 


Glycls 


PMP1 


28 


Extcl. Trnsp 


SER2 


28 


GlySer 


SER1 


28 


GlySer, Pyrldoxine 


RPE1 


27 


PentosePhosphate 


GCV3 


23 


Folate, GlySer 


LPD1 


23 


TCA cycle, GlySer, Glycls 


ACH1 


14 


Pyruvate 


ALD6 


12 


TyrTrpPhe 


FUM1 


11 


OxPhos 


A ATI 


3 


AlaAsp , ArgPro , TyrTrpPhe 


APA1 


2 


Cys , PurPyr 


IDP1 


2 


TCA cycle 


ALT2 


2 


AlaAsp 


GLY1 


2 


ThrLys 


PSD1 


2 


Plip 


MET22 


2 


Cysteine , 


PROl 


1 


ArgPro 


DCD1 


1 


Nuc . Salv . PurPyr 


ADOl 


1 


Nuc . Salv . 


PR02 


1 


ArgPro 



Table S5. Pleiotropy is present in catabolic genes and absent in biosynthetic ones. 
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Full name 

Alternate Carbon Metabolism 
Transport Mitochondrial 

Tyrosine Tryptophan and Phenylalanine Metabolism 

Arginine and Proline Metabolism 

Threonine and Lysine Metabolism 

Transport Extracellular 

Other Amino Acid Metabolism 

Glutamate metabolism 

Valine Leucine and Isoleucine Metabolism 

Pyruvate Metabolism 

Citric Acid Cycle 

Nucleotide Salvage Pathway 

Folate Metabolism 

Phospholipid Biosynthesis 

Purine and Pyrimidine Biosynthesis 

Cysteine Metabolism 

Glycine and Serine Metabolism 

Methionine Metabolism 

Alanine and Aspartate Metabolism 

Glycerolipid Metabolism 

Arabinose Metabolism 

Oxidative Phosphorylation 

Glycolysis and Gluconeogenesis 

Pentose Phosphate Pathway 

Anaplerotic reactions 

Pyridoxine Metabolism 

Table S6. Notation of metabolic modules at 
used thoughout the supplement. 



Abreviated name 

Alt . CS 
Mit . Trnsp 
Tyr,Trp,Phe 
Arg, Pro 
Thr, Lys 
Extcl . Trnsp 
OAA 
Glu 

Val,Leu,Ile 
Pyr 

TCA cycle 

Nuc. Salv. 

Folate 

Plip 

Pur , Pyr 

Cys 

Gly,Ser 
Met 

Ala, Asp 
Glip 
Ara 
OxPhos 
Glycls 
Pent . Ph . 
Anapl . 
Pdx 

jciated to Fig. 2, main manuscript, that is also 
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Figure SI. Differential stability of interactions: influence of metabolic 
distribution, strength, and sign. Differential stability of WT interactions with respect to 
whether their constituents genes belong to the same metabolic module (intra-module vs. 
inter- module; left), interaction strength (strong vs. weak; center) and type (positive vs. 
negative; right). For each interaction, instability was computed as the number of backgrounds 
where it disappears, changes sign, strength or both. In this figure, bars represent averages 
(p- values shown on top were obtained using the wilcoxon test between each pair of groups) . 
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Figure S2. Genetic hubs in the WT iND750 metabolism of S. cerevisiae. We depict here the 
degrees of the top 25 most connected genes (this is the list of genetic hubs considered in the 
manuscript). Colors indicate the proportion of negative and positive interactions. 
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Figure S3. Schematical representation of the metabolic changes that occur upon MIR1 
deletion. Colors indicate whether a given flux has increased, decreased, shut off or on, or 
changed direction, as shown in the legend. 
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Figure S4. Rewiring of genetic interactions in MIR1 background (AMIRl). Top: disrupted 
negative (yellow) and positive (blue) interactions. Node color represents the metabolic effect of 
the deletion on the reaction(s) associated to the node (see legend). Almost all the lost 
interactions involve either the background gene or genes that interact positively with it. 
Bottom-left: new negative interactions. Bottom-right: new positive interactions. 
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GR=0.42 



Figure S5. A) Schematic summary of metabolic and genetic changes after a background 
deletion. The cartoon is based on the AMIR1 background, but aiming to provide a description 
of the general principles of genetic rewiring, and how it reflects metabolic changes. A) Upon 
deletion of a background gene (BG), a compensatory mechanism is activated (Cl). This 
compensation usually originate side effects (in the case of MIR1, extrusion of malate from 
mitochondria) unless it is fully equivalent (e.g., duplicated genes). Additional mechanisms are 
thus activated to correct these secondary imbalances (Cl, C2, C3, and C4, dashed lines). B) 
WT genetic network, with nodes associated to flux-carrying reactions colored green, and single 
fitness contribution as node size. BG interacts negatively (yellow) with compensatory 
mechanisms (Cl, C2, C3, C4). Positive interactions (blue) occur either with the targets of its 
function (e.g., with respiration, R, as it is the ultimate target of phosphate import to 
mitochondria) or with genes that stop functioning when BG is compromised (generally 
displaying lower fitness effect thang BG; SP for strong positive partner). C) In the mutated 
network, the target mechanism (respiration) interacts positively with those genes that 
compensate for the BG deletion (Cl, C2). These compensatory mechanisms can interact among 
them positively or negatively, depending on the underlying functional relationship. Some of 
them can have in turn their own compensating mechanisms in the new background (e.g., C'l 
buffers Cl) also producing new negative interactions. 





Figure S7. Rewiring of hubs GCV3 (A) and SER2 (B) in all backgrounds. Figures equivalent 
to Fig. 4A, main text. 





Figure S9. Rewiring of hubs C0X1 (A) and PGI1 (B) in all backgrounds. Figures equivalent 
to Fig. 4A, main text. 
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Figure S10. Rewiring of hubs TPI1 (A), FBAl (B) and IDP1 (C) in all backgrounds. 
Figures equivalent to Fig. 4A, main text. 




Figure Sll. Rewiring of hubs ALT2 (A) and AAT1 (B) in all backgrounds. Figures 
equivalent to Fig. 4A, main text. 
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Figure S13. Rewiring of hubs RPE1 (A) and ALD6 (B) in all backgrounds. Figures 
equivalent to Fig. 4A, main text. 
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Figure S14. Rewiring of hubs GLY1 (A), SDH2 (B), IDP1(C) and IRC7 (D) in all 
backgrounds. Figures equivalent to Fig. 4A, main text. 
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Figure S15. A) Distribution of the connectivity (degree) exhibited by each hub in each 
background. Red dots show the WT degree, thick gray lines illustrate percentile 10 to 90, and 
dashed lines represent the range from the minimum (nonzero) to the maximum degree values 
observed. B) Number of gained and (number of) conserved genetic interactions are negatively 
correlated in hubs (each dot represents the Spearman's correlation coefficient between lost WT 
interactions and newly gained ones; red color indicates P <0.001 significance). C) We classified 
the type of rewiring as {horizontal: 23 hubs, vertical: 37 backgrounds where they become 
unstable): "Collapse", if the hub lost all its interactions; "Loss", if it lost only part; "Partial 
recovery" if part of the missing genetic interactions were substituted by new ones; "Dynamic 
rewiring" if all of the lost interactions were substituted by new ones, so that the new degree is 
at least equal to that of the WT; "Gain" if the hub conserves the same partners as in the WT 
but gained additional ones. Moreover, "deleted" means that the node is the same as the 
background, and "conserved" means that the interactions of the node remain the same. 
Background names are colored black if they are as well hubs, red if they are not hubs in the 
WT but become a hub in at least one background, and gray if neither of these definitions apply. 
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Figure S16. Distribution of the number of nodes and interactions, by interaction type, in the 
200 genetic networks resulting from neutral deletion accumulation trajectories. Red triangles 
indicate WT values. 



27 



A 



o 



CM 



B 



SP 



WP ■ 




E 
o 



SL 



N 



Disappear Appear 



WN WP 

To 



Transitions 

20 




Figure S17. A) Average number of negative (N) and positive (P) interactions observed to 
disappear and appear in genetic networks in response to neutral backgrounds. B) Average 
number of transitions between interaction types observed when comparing these networks with 
the WT one. A stronger tendency to rewiring can be observed in negative interactions overall. 
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Figure S18. Differential rewiring in LG and HG backgrounds. We show (left to right): 
average number of nodes whose degree is conserved (all partners equal than in the WT), 
number of those that conserved degree but changed their partners, number of nodes that lost 
interactions, and number of nodes that gained interactions. 



Network size (# Nodes) Network size (# Interactions) Mean degrees 

P = 1e-8 P = 2.7e-7 P = 0.01 




30 40 50 60 70 80 90 100 200 300 400 12 3 4 5 



Number of nodes Number of edges # Interactions / # Nodi 



# SL Nodes 
P = 0.39 




30 40 50 60 70 
Number of nodes 



# SL Interactions 
P = 0.84 




Mean SL degree 
P = 0.89 




0.2 0.4 0.S 0.8 1.0 1.2 

# Interactions / # Nodes 




#WP Nodes 
P = 2.6e-7 



# WN Interactions 
P = 2.5e-10 




Mean WN degree 
P = 0.18 





#WP Interactions 
P = 3.4e-8 



Number of edges 



Mean WP degree 
P = 0.001 




1 5 2 0 2.5 3.0 : 

# Interactions / # Nodes 




D 15 20 25 30 
Number of nodes 



#SP Interaction 
P = 3.37e-7 




Mean SP degree 
P = 0.001 




3 1.5 2 0 2.5 

# Interactions / # Nodes 



Figure S19. Top: distribution of number of interactions, number of nodes, and 
interactions/node for LG (red) and HG (blue) networks (neutral backgrounds). Bottom: 
before but restricted to the subnetworks constituted by synthetic lethal (SL, left), weak 
negative (WN, middle-left), weak positive (WP, middle-right), and strong positive (SP, ri 
interactions. 
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Figure S20. Frequency of LG and HG neutral trajectories in which a given WT interaction 
(dot) rewires (i.e., it dissappears, changes strength or sign). A) For each interaction, a p- value 
was obtained using Fisher's exact test (comparing the number of instabilities observed in HG 
and LG) and correction for multiple testing using the Benjamini-Hochberg-Yekutieli 
procedure [2]; those cases with p < 0.001 were colored in orange. B) Fraction of weak and 
inter- module interactions in the "LG-enriched" and "rest" subgroups. 





Figure S21. A) Distribution of the number of new interactions by type in LG (L) and HG 
(H) genotypes. For reference, we plotted the distribution of interactions in the WT network as 
a gray histogram line. B) For all unique interactions that appear in any of these genotypes, we 
plot its frequency in the HG and LG classes. We colored in red those interactions significantly 
enriched in LG networks (individual Fisher tests, p- values corrected for multiple testing using 
the Benjamini-Hochberg-Yekutieli procedure [2]). C) Interactions shown in red in B are 
represented here as a network (edge thickness is proportional to the frequency of each 
interaction). D) Transition map between types of interactions in LG, and E) in HG networks. 
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Figure S22. A) For each node in the WT network, we computed the frequency it disappears 
from the network in LG and HG genotypes (in red those nodes that are lost more frequently in 
LG). B) The same was done measuring the frequency at which genes become essential in LG 
and HG genotypes (in red those genes that become more frequently essential in LG). 
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Figure S23. Correlation between genetic rewiring in different backgrounds and their A) 
degree, B) total flux through the associated reactions, C) pleiotropy, and D) single-gene fitness 
contribution. 
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Figure S24. We selected from the high-confidence genetic interaction dataset in [3] those 
interactions between genes present in the iND750 model of S. cerevisiae. We separated the 
dataset into strong positive (upper quartile among positive in strength), strong negative (idem 
among negative) and weak (the rest). For functional classification, we conservatively 
considered a gene as "catabolic" if it belongs to either glycolysis, TCA cycle or oxidative 
phosphorylation; the rest were considered "biosynthetic" . According to this classification, an 
interaction can be catabolic, biosynthetic, or mixed (when one gene is catabolic and the other 
biosynthetic). In each strength category, we computed the percentage belonging to each 
functional class (vertical lines). To assess if any class is particularly enriched or depleted, we 
performed a bootstrapping analysis by randomizing 10000 times the assignment between 
strength and functional category of interactions. The distribution of percentages is shown in 
histograms, with the corresponding p-values. Plot boxes were colored red if significant 
depletion, and green if significant enrichment. 
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Figure S25. We computed, for all interactions appearing in at least one treatment, but not in 
the untreated network (i.e., treatment-specific interactions), the proportion of each 
strength/sign category. Treatment specific interactions are enriched in weak (especially weak 
positive) and depleted in strong (especially, in strong negative). 
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Figure S26. Genetic interactions between three genes representative of pentose phosphate 
pathway (ZWF1), glycolysis (PGI1) and TCA cycle (LPD1) in genotypes resulting from HG 
and LG trajectories. Each row represents the interactions associated to one genotype; in both 
cases the bottom row shows the epistasis in the WT for reference. It can be observed that the 
WT structure is conserved in a significant proportion of HG genotypes. In contrast, LG 
genotypes show a) less conservation of the WT pattern, b) generally more epistasis than in the 
WT, and c) a predominance of negative interactions, all of them patterns also significant at 
the statistical level (main text). 
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Figure S27. Synthetic lethal cluster 1. 

Synthetic lethal interactions in Saccharomyces cerevisiae iND750 model 

We inspected in detail the mechanistic underpinnings of synthetic lethal (SL) interactions in the 
WT iND750 metabolic model of S. cerevisiae. In contrast to the rest of interactions (which form 
a densely connected component), SL ones appear dispersed in several smaller network clusters 
that are mostly associated to biosynthetic functions. 

Synthetic lethal cluster 1 

We start with the biggest cluster shown in Fig. S27. The interactions that compose this cluster 
are specified in the following table, 
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Gene Syst. name Reactions coded 



ALT2 YDR111C Alanine transaminase 

a-ketoglutarate + Ala ^ Glu + pyruvate 

BNA2 YJR078W Tryptophan 2 , 3-dioxygenase 

Trp + O2 ^ N-f ormyl-kynurenine 

BNA5 YLR231C Kynureninase 

kynurenine + H20 ^ anthranilate + Ala 

SER2 YGR208W Phosphoserine phosphatase 

phosphoserine + H20 ^ Ser + Pi 

SER1 Y0R184W3 Phosphoserine aminotransferase 

phosphoserine + a-ketoglutarate ^ 3-phosphohydroxypyruvate + Glu 

AGX1 YFL030W Alanine : glyoxylate aminotransferase 

Ala + glyoxylate ^ pyruvate + glycine 

ICL1 YER065C Isocitrate lyase 

isocitrate ^ glyoxylate + succinate 

CTP1 YBR291C Mitochondrial inner membrane citrate transporter 

citrate [c] + malate[m] ^ citrate [m] + malate[c] 

AAT1 YKL106W Mitochondrial aspartate aminotransferase 

a-ketoglutarate + Asp ^ oxaloacetate + Glu 

Mitoch. erythrohydroxyglutamate a-ketoglutarate aminotransferase 
Mitoch. tyrosine aminotransferase 

AAT2 YLR027C Aspartate aminotransferase 

a-ketoglutarate + Asp ^ oxaloacetate + Glu 
aspartate aminotranf erase peroxisomal 

Erythrohydroxyglutamate a-ketoglutarate aminotransferase 

Perox. erythrohydroxyglutamate a-ketoglutarate aminotransferase 

tyrosine aminotranf erase 

a-ketoglutarate + Asp ^ oxaloacetate + Glu 
Perox. tyrosine aminotranf erase 

IDP2 YLR174W Isocitrate dehydrogenase (NADP+) 

isocitrate + NADP+ ^ a-ketoglutarate + C02 + NADPH + H+ 



Here, the interactions between ALT2 and both BNA2 and BNA5 reflect alternatives for alanine 
biosynthesis. In WT conditions, alanine is synthesized by ALT2 from pyruvate and glutamate. 
In absence of ALT2, Alanine starts to be synthesized through the Trp degradation pathway via 
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Figure S28. Alternative pathways for synthesis of Ala, Ser, Gly and Trp. Gray arrows 
indicate fluxes in the WT metabolism. Blue arrows indicate fluxes that appear when gene 
names colored in blue are deleted; analogously, green arrows are fluxes that activate upon 
deletion of the gene colored in green. 

kynureine, for which both BNA2 and BNA5 are needed. In this path, Ala is synthesized with 
anthranilate as a byproduct of the kynureinase reaction. Anthranilate returns again to Trp 
production; and Trp reenters the cycle to continue producing Ala. The excess of Trp needed to 
compensate its outflow to Ala is supplied by higher flux in Trp synthesis from chorismate, see 
pathways in the following web links: 

tryptophan degradation via kynurenine 

tryptophan biosynthesis 

superpathway of phenylalanine, tyrosine and tryptophan biosynthesis 

Moreover, the interactions between SER1 and SER2, and CTPl, ICLl and AGX1 reflect alterna- 
tives for Ser and Gly production (as the latter is produced from Ser). In the WT, Ser is produced 
from glycolytic 3-phosphoglycerate by SER1 and SER2. In absence of any of these genes, Gly 
starts being synthesized from glyoxilate (for which CTPl, ICLl and AGX1 are required) and Ala. 
These metabolic changes are delineated in Fig. S28. 

Fig. S28 further shows that, if both ALT2 and one of SER2 and SER1 are lacking, Ala is needed 
to obtain Gly, Ser and Trp, but Ser and Trp are needed for Ala, originating the SL interactions 
among these three genes. Indeed, we were able to recover growth in these double mutants when 
supplemented with one of Ala, Gly or Ser. Supply of Trp on the contrary did not recover growth 
because Ser is synthesized from Trp but is in turn needed to refill the Trp degradation cycle 
through kynureine. 
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In this first SL cluster (Fig. S27), we identify as well interactions between AAT1, AAT2, CTP1 
and IDP2. The Aspartate aminotransferases AAT2 (cytosolic) and AAT1 (mitochondrial) consti- 
tute in this model the only alternatives for Asp biosynthesis. But some additional metabolic 
readjustments occur upon removal of AAT1. Firstly, Glu produced in mitochondria from a- 
ketoglutarate has to be imported from the cytosol in order to synthesize 2-oxoadipate, in the 
Lys synthesis pathway. As the system starts to consume more Glu than synthesized, this be- 
comes restored through export of mitochondrial citrate to cytosol (by CTPl), and transformation 
to a-ketoglutarate (by IDP2; a-ketoglutarate is converted in turn to Glu). This underlies the 
SL interactions of A ATI with CTPl and IDP2. Indeed, in the CTP1-AAT1 or IDP2-AAT1 double 
deletion, growth can be recovered by adding either Glu or a-ketoglutarate. Moreover, if the GDH3 
gene (a-ketoglutarate to Glu) is removed, a-ketoglutarate stops being able to recover growth in 
this double deletion, showing that its supply is used for replacing Glu. 
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Figure S29. Synthetic lethal cluster 2. 



Synthetic lethal cluster 2 

In the second cluster (Fig. S29), SL interactions with GLY1 denote different alternatives for Thr 
biosynthesis. In the WT, Thr is synthesized from Gly and acetaldehyde by GLY1. Alternatively, 
Thr can be also obtained from homoserine through THR1 and THR4. In turn, H0M2, H0M3 and HDM6 
are needed for biosynthesis of homoserine from Asp, and this is the rationale of the interactions 
of these genes with GLY1. 

Interactions in the left of Fig. S29 are linked to homocysteine biosynthesis (needed for Met 
and phosphatidilcholine) . In the pathway used in the WT, Thr is transformed to 2-oxobutanoate 
(in the Thr deaminase reaction, coded by two redundant genes); 2-oxobutanoate is then trans- 
formed into cystathionine by YML082W (and succinylhomoserine lyase, a gene-orphan reaction in 
the model). Finally, cystathionine is transformed into homocysteine by IRC7. In the alternative 
pathway, homocysteine is synthesized from Asp via homoserine in a pathway involving H0M3, 
H0M2, H0M6, MET2 and MET17. These being the two alternatives for Met biosynthesis; thus, genes 
in one pathway interact negatively with genes in the other. 
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Gene 



Syst. name Reactions coded 



GLY1 YEL046C Threonine aldolase 

Thr ^ Gly + acetaldehyde 

THR4 YCR053W Threonine synthase 

phosphohomoserine + H20 ^ Thr + Pi 

THR1 YHR025W Homoserine kinase 

homoserine + ATP ^ phosphohomoserine + ADP + 2H+ 

H0M2 YDR158W Aspartic beta semi-aldehyde dehydrogenase 

Aspartate semialdehyde + NADP+ + Pi ^ phosphoaspartate + H+ + NADPH 

H0M3 YER052C homoserine kinase activity 

homoserine + ATP ^ phosphohomoserine + ADP + 2 H(+) 

H0M6 YJR139C homoserine dehydrogenase 

homoserine + NADP+ ^ aspartate semialdehyde + NADPH + H+ 

IRC7 YFR055W cystathionine beta-lyase 

cystathionine + H20 ^ homocysteine + NH3 + pyruvate 

succinylhomoserine + H20 ^ 2-oxobutanoate + succinate + NH4 + H+ 

MET2 YNL277W homoserine Q-transacetylase 

homoserine + acetyl-CoA ^ Acetylhomoserine + CoA 

MET 17 YLR303W Q-acetylhomoserine (thiol) -lyase 

Acetylhomoserine + H2S ^ homocysteine + acetate 



YML082W YML082W 



An additional question is why, while MET2 interacts with both YML082W and IRC7, MET 17 
interacts only with IRC7. The explanation is that cystathionine can be alternatively synthesized 
from acetylhomoserine and Cys by STR2 (which is a slightly more inefficient mechanism). See 
also the web links below 

superpathway of threonine and methionine biosynthesis 

threonine degradation 

superpathway of methionine biosynthesis 
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Figure S30. Synthetic lethal cluster 3. 

Synthetic lethal cluster 3 

Gene Syst. name Reactions coded 

PROl YDR300C Gamma-glutamyl kinase 

Glu + ATP ^ glutamyl-5-phosphate + ADP + H+ 

CAR2 YLR438W L-ornithine transaminase 

ornithine + a 2-oxoacetate ^ glutamate semialdehyde + Glu 

PRQ2 Y0R323C Gamma-glutamyl phosphate reductase 

glutamate-5-semialdehide + NADP+ + Pi ^ glutamyl-5-phosphate + H+ + NADPH 

In WT conditions, Pro is synthesized from Glu in a linear pathway that sequentially involves 
PROl and PRD2 whose result is glutamate-5-semialdehyde. When one of these genes is mutated, 
an alternative for Proline biosynthesis is to use CAR2 in the Arg degradation pathway. See also 
web links below 

arginine degradation 

proline biosynthesis 
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Figure S31. Synthetic lethal cluster 4. 

Synthetic lethal cluster 4 

Gene Syst. name Reactions coded 

PSD2 YGR170W phosphatidylserine decarboxylase, golgi and vacuole 

phosphatidylserine + H+ ^ phosphatidylethanolamine + C02 

PSD1 YNL169C Phosphatidylserine decarboxylase, mitochondria 

phosphatidylserine + H+ ^ phosphatidylethanolamine + C02 

ATH1 YPR026W acid trehalase (vacuole) 

trehalose + H2Q ^ 2 glucose 

In WT conditions, phosphatidilserine is synthesized in mitochondria by PSD1. An alternative 
mechanism is its synthesis in vacuole by PSD2. For this, however, trehalase has to be imported 
to vacuole and degradated to glucose, supplying the vacuolar proton needed for PSD2. Therefore, 
both enzymes are SL with PSD1. 
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Figure S32. Synthetic lethal cluster 5. 

Synthetic lethal cluster 5 

Gene Syst. name Reactions coded 

TPIl YDR050C triose-phosphate isomerase 

glyceraldehyde-3-phosphate ^ glycerone phosphate 

PFKl YGR240C phosphof ructokinase 

ATP + f ructose-6-phosphate ^ ADP + f ructose-1 ,6-bisphosphate 

FBPl YKL060C Fructose 1 ,6-bisphosphate aldolase 

f ructose-1 ,6-bisphosphate ^ DHAP + glyceraldehyde-3-phosphate . 

Apparently, these enzymes should not interact negatively, as they perform their function se- 
quentially in glycolysis. However, they present an additional function necessary for growth, the 
generation of DHAP in a small, but essential, quantity that is required for phospholipid biosyn- 
thesis. DHAP can only be generated by TPIl from glyceraldehide-3-phosphate, or by FBPl 
from fructose- 1,6-bisphosphate which is in turn produced by PFKl. Thus, although most of the 
DHAP from the fructose bisphosphatase reaction is transformed to glyceraldehyde-3-phosphate 
for glycolytic progression, a very small quantity of DHAP is essential for biosynthetic causes. 
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Figure S33. Synthetic lethal cluster 6. 

Synthetic lethal cluster 6 

Gene Syst. name Reactions coded 

ALD6 YPL061W aldehyde dehydrogenase 

acetaldehyde + NADP+ + H2Q ^ acetate + NADPH + H+ 

ALD4 Y0R374W Mitochondrial aldehyde dehydrogenase 

acetaldehyde + NADP+ + H2Q ^ acetate + NADPH + H+ 

ADH3 YMR083W Mitochondrial alcohol dehydrogenase 

ethanol + NAD+ ^ acetaldehyde + NADH + H+ 

These interactions reflect the two alternatives (cytosolic and mitochondrial) to obtain acetate, 
which is in turn needed to produce acetyl-CoA. 

Rest of synthetic lethal interactions (isolated pairs) 
MET22 - PMP1 

Gene Syst. name Reactions coded 

PMP1 YCR024C-A ATPase cytosolic 

ATP + H20 -> ADP + Pi + H+ 

MET22 Y0L064C 3,5-bisphosphate nucleotidase 

adenosine-3' , 5 ' -bisphosphate + H20 ^ adenosine 5'-phosphate + Pi 



PMP1 codes in the model for ATP hydrolysis to produce ADP, important for adjusting the ATP 
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f± ADP homeostasis. As an alternative, adenosine-5-phopsphosulfate is transformed to AMP 
by MET22, that combines later with ATP to produce 2 ADP by adenilate kinase (a flux that also 
increases considerably in APMPl). This alternative pathway to ADP homeostasis maintenance 
consumes 2 extra ATP molecules, which is the explanation of the fitness decrease on APMPl. 

APA1 - MET3 

Two ways for the first step of S04 fixation to produce adenosine-5-phopsphosulfate. Through 
APA1 -used in WT conditions- it is fixed over ADP; through MET3, it is fixed over ATP and 
produces pirophosphate that is an energetically more costly alternative. 

AC02 - AC01 

Mitochondrial aconitase redundancy. Additionally, AC01 also performs cytosolic aconitase, 
which is the reason of keeping both genes in the dataset (see Materials and Methods, main text). 

ACH1 - ACS2 

Two different ways to synthesize acetyl-CoA from acetate. One of them involves ATP hy- 
drolysis to AMP and PPi (ACS2), and the other does not (ACHl) and is the alternative used by 
the WT. The energetic difference is responsible for a lower fitness of A ACHl 

LEU4 - LEU9 

2 isopropylmalate synthase redundancy. Additionally, LEU4 can also perform this activity in 
mitochondria, this difference being responsible for keeping both in the dataset. 

TPS3 - TSL1 

These two genes both code for two reactions, trehalose phosphatase and alpha-trehalose 
phosphate synthase. They are maintained in the dataset because of a different role in these 
reactions: while TSL1 is an essential sububnit, TPS3 is interchangeable (by TPS2, which is not 
kept in the dataset). 
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AR03 - AR04 

Both genes code for 3-deoxyarabinoheptulosonate-7-phosphate synthetase. However, ARD3 
codes additionally for the mitochondrial version of the reaction. 

Synthetic lethal interactions: concluding remarks. 

Synthetic lethal interactions were analyzed in this supplementary section in order to better 
understand their mechanistic underpinnings. Generally, most SL interactions arise from either 
paralell pathways or redundant genes. Some more intrincate architectures (such as the one 
associated to synthesis of Ala, Trp, Ser and Gly, see Figs. S27- S28) were also found. In that 
particular case, a set of inter-conversions between several essential biomass constituents led to 
a relatively big cluster of SL interactions. 

It is also important to emphasize that SL interactions are systematically associated to biosyn- 
thetic functions. This is not observed only for genes that clearly belong to biosynthetic pathways. 
For instance, although TPI1, PFK1 and FBA1 (Fig. S32) are clearly catabolic genes, the SL in- 
teractions between them reflect an important biosynthetic function that is linked to providing 
DHAP for phospholipid biosynthesis. Moreover, these genes also display other types of inter- 
actions (e.g., weak negative or positive) which are clearly related to their catabolic activity. 
This further stresses the functional segregation of SL and the rest of epistatic interactions in 
biosynthetic and catabolic functions, respectively. 
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